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1.  SUMMARY 


In  gallium  nitride  (GaN)-based  high  electron  mobility  transistor  (HEMT)  devices,  the  heat  flux 
could  reach  a  few  kW/cm2  and  the  heat  dissipation  has  become  a  critical  issue.  Due  to  the 
nanometer  length  scale  of  the  device  layers  and  the  high  density  of  interfaces,  non-diffusive  heat 
conduction  plays  a  critical  role  in  thermal  transport  of  GaN  devices  manifested  as  thermal 
conductivity  and  the  interfacial  thermal  resistance  of  the  device  layers.  Therefore,  understanding 
the  phonon  transport  physics  and  reducing  the  near-junction  thermal  resistance  are  of  critical 
need  for  practitioners.  This  work  uses  both  theoretical  and  experimental  method  to  investigate 
the  non-diffusive  phonon  transport.  From  density  functional  theory  (DFT)  based  Boltzman 
transport  equation  (BTE),  phonon  mean  free  paths  of  GaN,  silicon  carbide  (SiC),  and  diamond 
are  predicted  to  span  in  a  wide  range,  from  dozens  ofnanometers  to  several  microns.  Such  wide 
range  of  mean  free  path  (MFP)  indicates  that  the  non-diffusive  phonon  transport  would  even  be 
significant  at  room  temperature.  We  then  use  an  ultrafast  pump-probe  thermoreflectance 
measurement  system  to  characterize  the  thermal  conductivity  of  GaN  layers,  and  we  found  that 
the  measured  thermal  conductivity  is  related  to  the  length  scale  of  temperature  gradient  (heating 
frequency)  and  the  thickness  of  the  thin  film  heaters.  Both  our  theoretical  and  experimental 
work  suggest  that  the  thermal  dissipation  would  be  much  less  efficient  than  the  bulk  case,  due  to 
the  high  frequency  transient  heat  transport.  This  research  suggests  that  the  size  effect  of  thermal 
conductivity  in  GaN  devices  is  collectively  determined  by  several  characteristic  lengths,  and  the 
previous  knowledge  of  size  effect  is  too  simple  to  understand  heat  dissipation  in  practical  HEMT 
devices. 
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2.  INTRODUCTION 


GaN-based  HEMTs  are  promising  for  next  generation  radar  applications  due  to  their  wide  band 
gap,  high  electron  velocity,  and  high  breakdown  field.  However,  due  to  the  high  power  density 
in  GaN  HEMT  devices,  localized  self-heating  with  a  heat  flux  approaching  tens  ofkW/cm2 
can  be  a  critical  issue.  Figure  1  shows  a  state-of-the-art  GaN  HEMT  device  geometry,  where 
AlGaN/GaN  layers  are  heteroepitaxially  grown  on  non-GaN  substrates,  such  as  high  thermal 
conductivity  SiC  (~400  W  m  1  K  '  )  and  chemical-vapor-deposited  (CVD)  diamond 
(1000-2000  W  m_1  K1)  substrates,  due  to  the  limited  commercial  availability  of  bulk  GaN 
substrates.  Due  to  the  lattice  mismatch  between  GaN  and  SiC/diamond,  there  is  a 
transition/adhesion  layer  from  a  couple  of  microns  to  tens  of  nanometers  (latest  generations, 
usually  proprietary)  between  GaN  and  high  thermal  conductivity  SiC/diamond  substrates.  These 
transition/adhesion  layers  are  usually  amorphous  materials  with  defects,  such  as  AlGaN,  AIN,  or 
SiN,  which  can  be  the  most  resistive  thermal  pathway  in  the  device  and  thus  offsetting  the 
advantages  of  using  high  thermal  conductivity  substrates.  A  key  research  task  in  high  power  GaN 
HEMTs  is  to  reduce  the  “near-junction”  thermal  resistance,  as  identified  by  the  Defense 
Advanced  Research  Projects  Agency  (DARPA)  Microsystems  Technology  Office  (MTO). 


Al 


plasma  oxide  (1-2  nm) 


In  AIN  (7.4  nm) 


AIN  (1.0  nm) 


GaN  Channel  (26  nm) 


InGaN  back-barrier  (3.3  nm) 


GaN  Buffer  (l.a  um) 


SiC  substrate 


Figure  1:  Schematic  of  GaN  HEMT  Device  Geometry 

To  study  the  near-junction  thermal  transport  and  to  reduce  the  near-junction  thermal  resistance  of 
GaN  devices,  new  strategies  must  be  created  to  quantify  contributions  of  each  phonon  modes  to 
the  thermal  conductivity  of  the  material  stack  in  GaN-on-SiC  and  GaN-on-diamond  devices. 
Measurement  of  phonon  mean  free  path  distributions  in  the  materials  at  different  temperature 
appears  to  be  the  most  relevant  yet  most  challenging  task  that  can  potentially  bridge  the  rapid¬ 
growing  fundamental  understanding  of  phonon  transport  with  practitioners’  interest  in  device 
scaling. 

As  a  statistical  concept,  the  phonon  MFP  measures  the  average  travelling  distance  between 
two  consecutive  phonon  scattering  events.  Phonon  MFP  depends  strongly  upon  the  phonon  modes 
and  scattering  details.  The  phonon  scattering  events  include  the  3 -phonon  and  4-phonon  scattering 
due  to  the  anharmonic  nature  of  the  interatomic  potential,  phonon-defects,  phonon- 
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boundary/interface  scattering.  Contributions  of  phonons  with  different  MFPs  to  heat  transfer  can 
be  examined  through  the  lattice  thermal  conductivity  k  along  a  direction  predicted  by  the  kinetic 
transport  theory1 : 


b  = 


-Z 


C  va  Aa 

^qs  ‘/qsilgjs 


<1* 


(1) 


where  q,s  represents  the  5-th  phonon  mode  at  q  in  the  first  Brillouin  zone,  Cqs  is  the  mode- 
specific  heat  capacity,  vqs  is  the  group  velocity,  and  the  MFP  Aqs  for  each  phonon  mode  is  the 
product  of  spectral  phonon  group  velocity  and  relaxation  time  zqs,  i.e.,  Aqs  =  vqsTqs.  For  a 
perfect  bulk  crystal,  phonon  MFP  spans  several  orders  of  magnitude,  from  a  few  nanometers  to 
tens  of  microns.  The  thermal  conductivity  can  therefore  be  written  as  an  integration  over  the 
entire  phonon  MFP: 


k  = 


f 


<p(A)dA 


(2) 


If  the  phonon  mean  free  path  is  comparable  or  even  smaller  than  the  characteristic  length  of  the 
device  Lc,  the  thermal  conductivity  would  reduce  due  to  the  suppressed  phonon  MFP: 


kA~*  [  <p(A)dA 

Jo 


(3) 


where  /ca  is  the  apparent  thermal  conductivity  when  the  phonons  transport  is  non-diffusive. 

The  wide  range  of  phonon  spectrum  suggest  that  the  thermal  transport  in  HEMT  nano-devices 
would  significantly  deviate  from  the  diffusive  transport  in  bulk  limit  due  to  the  boundary 
scattering  with  the  interface.  Besides  the  boundary  length  of  the  GaN  layer,  there  are  also  other 
length  scales  significantly  affecting  the  thermal  transport.  For  example,  the  length  scale  of  the 
nano-sized  heater2'3  and  the  length  scale  of  temperature  field4.  In  HEMT  devices  working  in 
GHz  range,  the  heating  and  thermal  dissipation  would  be  highly  transient.  We  can  therefore 
define  the  penetration  depth  as  length  scale  of  the  temperature  field: 


dp  yj  kjnCf  ^ 

where  /  is  the  heating  frequency.  However,  it  is  still  not  understood  how  different  length  scales 
simultaneously  determine  the  phonon  transport  regime.  Therefore,  in  this  project,  1)  we  first  use 
first  principles-based  calculations  to  predict  the  phonon  mean  free  path  spectrum  and  2)  we 
discovered  frequency-filtering  effect,  a  novel  phenomenon  that  the  thickness  of  the  metal  film 
and  heating  frequency  simultaneously  affect  the  measured  apparent  thermal  conductivity. 
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3.  FIRST  PRINCIPLES  CALCULATION  OF  PHONON  MFP  IN  GaN,  SiC 

AND  DIAMOND 


First  principle -based  methods  for  calculating  the  lattice  thermal  conductivity  for  a  variety  of  bulk 
and  nanostructured  materials  have  been  developed  recently,5"8  where  great  agreement  with  the 
measured  thermal  conductivity  values  at  different  temperatures  has  been  demonstrated.  These 
calculations  can  differentiate  the  contributions  to  the  total  thermal  conductivity  of  phonons  with 
different  wavelengths  and  different  mean  free  paths.  Here,  we  calculate  the  thermal  conductivities 
and  phonon  mean  free  path  distributions  of  diamond,  w-GaN  and  4H-SiC  (see  Figure  2)  using  the 
iterative  solution  of  Peierls-Boltzmann  transport  equation. 


Diamond 


w-GaN 


4H-SiC 


Figure  2:  Crystal  Structure  of  Diamond,  w-GaN  and  4H-SiC 


3.1  Peierls-Boltzmann  Transport  Equation 

To  calculate  the  thermal  conductivity  along  a  (=  a,  b,  c)  direction,  a  small  temperature  gradient 
was  applied  along  a  direction,  — ; ,  to  perturb  the  phonon  population  from  the  equilibrium 
Bose-Einstein  distribution,  rtq5  ,  to  the  non-equilibrium  one,  nqs.  When  the  steady  state  was 
achieved,  the  heat  flux  can  be  expressed  as  the  summation  of  the  contributions  from  all  phonon 
modes  through: 


hcoqsv%snqsdq 

(5) 

where  h  is  the  Planck  constant,  o)qs  is  the  frequency  of  phonon  mode  qs.  The  phonon  frequency 
o)qS  and  Vqs  can  be  obtained  from  the  phonon  dispersion  of  the  crystal,  while  nqs  is  determined 
by  the  Peierls-Boltzmann  Transport  Equation  (PBTE).  After  Ja  was  calculated  with  the 
contributions  from  all  phonon  modes,  the  macroscopic  thermal  conductivity  was  then  derived 
using  the  Fourier’s  law  of  heat  conduction,  Kaa  =  —]“/(—)  In  order  to  evaluate  the  heat  flux 
driven  by  the  applied  small  temperature  gradient,  non-equilibrium  phonon  distribution  of  each 
phonon  mode  in  Eq.  (5)  needs  to  be  determined. 
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The  phonon  dispersion  relation  is  related  to  the  second-order  harmonic  force  constants  of  the 
crystal,  <fi.  With  the  harmonic  force  constants,  the  dynamical  matrix  D  can  then  be  solved  for 
phonon  dispersion.  The  elements  of  D  are  written  as: 


<?(i)  = 


JM,Mr 


;X<Cv 


(6) 


where  Rr  is  the  r-th  basis  atom  in  the  unit  cell  which  is  represented  by  the  translation  vectors  R, 
Mr  is  the  atomic  mass  of  the  r-th  basis  of  the  primitive  cell.  The  phonon  frequency,  a)qs,  is  the 
square  root  of  the  s-th  eigenvalue  of  the  dynamical  matrix  and  the  group  velocity,  VqS,,  is 
calculated  as  d(Oqs/dqa. 


To  determine  the  non-equilibrium  phonon  population  distribution  nqs,  it  is  convenient  to  express 
nq s  in  the  form  of  nqi  -I-  rtqi  (nqv  +  l)^  with  deviation  function  X^.  The  linearized 
PBTE,  which  describes  the  balance  between  phonon  diffusion  driven  by  the  applied  temperature 
gradient  and  phonon  scattering,  was  then  solved  for  the  non-equilibrium  phonon  distribution 
function,  or  its  equivalent  jqs.  In  this  work,  we  consider  the  three-phonon  scattering  and  phonon- 
isotope  scattering,  and  the  corresponding  PBTE  can  be  expressed  as5'6: 


v. 


a 


q* 


dr 


=  Z  Klv  (4-V-  -  <V  - 


1  ci's'  a"s" 

+  -W 4  q 

2  ^ 


+ 


Z\wiso,  , 

L  * 

qV 


(*Jv'  + 
(<V  "  4*)] 


(?) 


The  first  two  terms  on  the  right  hand  side  of  Eq.  (4)  are  related  to  three-phonon  scattering  rates, 
and  the  third  term  is  the  scattering  rates  due  to  phonon-isotope  scattering.  Here,  s,  ,  and 
W* s ,q  5  are  expressed  as:  qi'q  i 


W 


qV  =  2im°snqV  (njv,  +  l)  |l4(-qs, -qV, q'V')l: 
x<S(ajqs  +  —  ttyf5fj)6(q  +  q'  —  q ''  +  G) 


(8) 


W, 


qVq'V1 

qs 


=  2irnqs  (n°v  +  l)  (n°,v,  +  l)  \V3(-qs, q' s', q"s")|2 

x5(cuqs  -  tcj qv  -  ^yvO  5(q  -  q'  -  q"  +  g) 


(9) 


where  the  8  function  and  the  A  function  denotes  the  energy  conservation  condition  (Dqs±  oiqv  - 
uiq "s"  =  0  and  the  momentum  conservation  condition  q  ±q’-q”  +  G  =  0  for  the  three-phonon 
scattering  process.  The  +  and  -  signs  represent  the  annihilation  and  decay  processes,  respectively. 
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The  G  vector  represents  a  reciprocal  vector.  The  scattering  is  normal  process  when  G  =  0,  while 
it  is  Umklapp  process  when  G  ^  0.  The  three-phonon  scattering  matrix  V3  is  written  as: 


VgCq^qV^'V')  = 


KSN0(i)qs(i>qrsia)qits" 


ip 


<*PY 

Ot.rV.R  "t" 


x  exp(tq'  ■  R'  +  iq''  ■  R") 


t  ftV  R*fTif  afiy 


eta  yp  pT>>y 

tqs  fcqVV's" 


(10) 


where  e  is  the  eigenvector  of  the  dynamical  matrix. 
The  isotope  scattering  matrix  jy  is  written  as: 


WlSOn,  ,  = 

qi  .q  a 


TT 


("Jv  +  *)  £  A  1(4") 


-tct 

-qV 


*( 


OJ, 


qs 


y) 


(11) 


where  g%  is  the  mass  variance  parameter  for  the  basis  atom  r9.  The  mass  variance  parameter  is 
related  to  the  concentration,  a,  and  the  mass,  Mi,  of  the  z-th  constituent  isotopes, 


9  = 


(12) 


In  order  to  sample  the  reciprocal  space  for  phonon  scattering,  and  thermal  conductivity 
calculations,  Eq.  (2),  we  discrete  the  first  Brillouin  zone  into  N\xNixN3  T-centered  q-points.  The 
grind  points  are  located  at  <J  =  —  bt  -I - b7  H - b3  where  bi,  b2,  and  b3  are  reciprocal 

Ni  N2  N3 

primitive  lattice  vectors,  and  to,  n  and  /  are  integers.  The  delta  function  for  the  energy 
conservation  condition  in  Eq.  (5)  and  Eq.  (6)  is  evaluated  by  the  Gaussian  function  with  an 
adaptive  broadening  parameter  related  to  the  local  group  velocity.10  The  set  of  linear  equations 
Eq.  (7),  with  respect  to  ifqi- ,  can  be  self-consistently  solved  through  the  iterative  method.  Here 
we  employ  the  bi-conjugate  gradient  stabilized  method,11  a  variant  of  the  conjugate  gradient 
algorithm,  to  iteratively  solve  Eq.  (7).  After  i£qy  is  calculated,  the  thermal  conductivity  can  be 
expressed  as 


^ aa  ^  ft^qjs^qi-ftqifftqs  ”1"  l)^qi 

Is  (13) 

To  quantify  the  importance  that  the  temperature  gradient  collective  scattering  behavior  of  all  the 
phonon  modes  plays  on  anisotropic  phonon  transport,  we  rewrite  Eq.  (10)  into  the  form  of  Eq. 

(1)  to  defining  the  phonon  MFP  in  the  iterative  approach,  which  reads  as 
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(14) 


%  =  (ksT2/h)  ■  (F‘/uv) 


where  ks  is  the  Boltzmann  constant. 

3.2  Extracting  Interatomic  Force  Constants  from  First-Principles 

All  first-principles  calculations  are  performed  using  the  Quantum-Espresso  package12  with  norm- 
conserving  pseudopotential  and  the  local-density  approximation  functional.  The  kinetic-energy 
cut-off  for  the  plane-wave  basis  is  set  at  80  Ry.  12x12x12,  6x6x6  and  8x8x4  k-meshes  are  used 
to  sample  electronic  states  in  the  reciprocal  space  of  diamond,  w-GaN  and  4H-SiC,  respectively. 
The  choice  of  the  energy  cutoff  and  /c-mesh  ensures  that  the  energy  change  is  smaller  than 
1  meV/atom  when  refining  these  two  parameters.  The  lattice  structures  are  relaxed  using  the 
conjugated  gradient  method  to  minimize  the  energy  of  the  crystal,  and  optimized  lattice  constants 
are  summarized  in  Table  1,  which  are  close  to  the  values  from  experimental  measurements.13"15 


Table  1.  Lattice  Constants  of  Diamond,  w-GaN,  and  4H-SiC 


Materials 

a  (nm) 

c  (nm) 

DFT 

Exp. 

DFT 

Exp. 

Diamond 

0.3527 

0.356713 

w-GaN 

0.3160 

0.319015 

0.5153 

0.518915 

4H-SiC 

0.3068 

0.307314 

1.0045 

1.005314 

To  calculate  the  lattice  thermal  conductivity  and  phonon  mean  free  path,  both  the  second-order 
harmonic  and  third-order  anharmonic  interatomic  force  constants  of  each  material  are  required. 
We  use  the  density-functional  perturbation  theory  (DFPT)  approach16  to  calculate  the  dynamical 
matrices  of  the  phonon  modes  on  the  6x6x6  q-meshes  for  diamond,  w-GaN  and  4H-SiC.  With 
the  dynamical  matrices,  the  harmonic  force  constants  are  calculated  through  the  Fourier- 
transform  of  the  obtained  dynamical  matrices.  Both  the  DFPT  and  Fourier- transform  are 
implemented  in  Quantum-Espresso. 


For  the  third-order  anharmonic  force  constants,  the  supercell  scheme  with  a  size  of  4x4x4, 

3x3x3  and  3x3x1  primitive  cells  of  diamond,  w-GaN  and  4H-SiC  is  employed.  We  displace  two 
atoms  in  a  supercell  by  A  =  0.01  A  away  from  their  equilibrium  position  along  ±x,  ±y  and  ±z 
directions  simultaneously,  and  then  record  the  forces,  F,  of  all  atoms  in  the  supercell.  With  the 
recorded  force  information  on  all  atoms  in  the  supercell,  the  third-order  force  constants  are 
calculated  using  the  finite-difference  scheme,17  which  are  expressed  as 


KSrVyV'  “  “  A'  u*  V  -  A)  “  Fl»t»{ulL  ~  -  -A)  - 


Fr"t "(“**  “  _A'U*V  -  A)  +  -  _A'U«V  -  -A)]> 


(15) 


where  u  is  the  atom  displacement  from  equilibrium  position.  The  Lagrange  multipliers 
approach18  is  employed  to  enforce  translational  invariance  for  the  third-order  force  constants, 
without  which  the  long-wavelength  phonons  are  likely  to  be  more  severely  scattered  and  thus 
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underestimate  the  thermal  conductivity.  The  cutoffs  for  the  third-order  anharmonic  interaction 
are  chosen  to  the  third  nearest  neighbors. 

3.3  Results 

Figure  3  shows  the  calculated  phonon  dispersion  curves  for  diamond,  w-GaN  and  4H-SiC.  We 
also  plot  the  phonon  dispersion  measured  by  neutron  scattering  experiments.  The  good 
agreement  with  both  the  measurements19"21  validates  the  interatomic  force  constants  from  our 
first-principles  calculations. 

The  calculated  temperature-dependent  thermal  conductivity  of  silicon  and  sapphire  is  presented 
in  Figure  4,  as  well  as  experimental  data  for  both  silicon22  and  sapphire23.  Our  calculations  are 
close  to  these  previous  works. 

The  accumulative  thermal  conductivity  gives  the  contribution  from  the  phonons  with  MFP 
smaller  than  A  to  the  total  thermal  conductivity.  By  summing  up  the  thermal  conductivity  of 
phonon  modes  with  MFPs  in  that  range,  the  accumulative  thermal  conductivity  can  be 
determined. 

Figure  5  shows  the  differential  and  cumulative  thermal  conductivity  of  diamond,  w-GaN  and 
4H-SiC  from  first-principles  calculations.  It  is  found  that  the  phonons  with  MFP  larger  than  1  um 
in  all  materials  studied  contribute  to  more  than  50%  of  the  total  thermal  conductivity.  Since  the 
thickness  of  GaN  thin  film  is  around  800  nm  in  HEMT  devices,  the  long-MFP  phonons  should 
behave  quite  different  from  the  prediction  from  the  Fourier’s  law  and  we  have  to  consider  how 
the  ballistic  transport  affects  the  heat  conduction  in  thin  films  as  well  as  phonon  transmission 
across  the  interfaces. 
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Frequency  (cm'1)  Frequency  (cm'1) 


Figure  3:  Calculated  Phonon  Dispersion  Relation  of  (a)  Diamond,  (b)  w-GaN,  and 

(c)  4H-SiC 

The  blue  dots  in  (a-c)  represent  the  measured  phonon  dispersion  of  diamond19,  w-GaN21,  and 

4H-SiC.20 
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Figure  4:  Thermal  Conductivity  of  (a)  Diamond,  (b)  w-GaN,  and  (c)  4H-SiC,  as  a  Function 

of  Temperature 

The  dash  lines  in  (a)  represent  the  calculated  thermal  conductivity  of  diamond  using  the  first- 
principles  based  PBTE  method 22 ,  while  the  green 23  and  blue24  dots  are  the  measured  thermal 
conductivity  of  diamond.  The  green25  and  blue26  dots  in  (b)  show  the  experimentally  measured 
thermal  conductivity  for  w-GaN.  The  dots  in  (c)  represent  the  measured  thermal  conductivity  of 

4H-SiC.27 
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Figure  5:  Cumulative  Thermal  Conductivity  as  a  Function  of  Phonon  Mean  Free  Path 
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4.  FREQUENCY  FILTERING  EFFECT  IN  TIME-DOMAIN 

THERMOREFLECTANCE 


Size  effect  in  thermal  conductivity  has  been  intensively  studied  in  the  past  two  decades  to 
understand  the  nanoscale  thermal  transport.  For  example,  the  thickness  of  GaN  layer,28  size  of 
the  heater,2'3, 29  and  thermal  penetration  length  (modulation  frequency).30"31  However,  the  size 
effect  of  thermal  conductivity  as  a  function  a  single  length  scale  is  over  simplified  for  thermal 
transport  in  HEMT  devices,  because  more  than  one  characteristic  lengths  in  HEMT  devices 
together  determines  the  phonon  transport  physics.  To  understand  how  different  length  scales 
simultaneously  affect  the  apparent  thermal  conductivity,  we  simply  the  multi-layered  HEMT 
device  to  a  two-layer  system  with  A1  thin  film  deposited  on  GaN  substrate.  The  A1  thin  film 
serves  as  the  heater  by  absorbing  the  pump  laser  to  create  a  periodic  heating.  Such  simplification 
reduces  the  number  of  characteristic  lengths  to  only  two:  the  thermal  penetration  depth  and  the 
thickness  of  the  A1  thin  film.  We  then  use  the  time-domain  thermoreflectance  (TDTR)  technique 
to  optically  measure  the  thermal  conductivity  of  GaN  in  the  Al/GaN  sample.  By  varying  the 
heating  frequency  and  the  thickness  of  the  A1  thin  film,  the  measured  apparent  thermal 
conductivity  is  changed.  We  explain  this  novel  phenomenon  as  the  frequency  filtering  effect  of 
the  A1  thin  film.  This  work  is  the  first  observation  on  how  different  feature  size  simultaneously 
determine  the  phonon  transport  regime. 

4.1  Implementation  of  the  TDTR  Measurement 

Time-domain  thermoreflectance  measurement  is  an  optical  experiment  to  characterize  thermal 
properties  of  materials.  Figure  6  shows  the  experimental  system  implemented  in  our  lab.  A 
Ti-sapphire  oscillator  is  used  as  the  light  source  with  wavelength  centered  at  around  800  nm  and 
a  repetition  rate  of  80  MHz.  The  laser  output  is  split  into  pump  beam  for  heating  and  probe  beam 
for  sensing.  Before  focused  onto  the  sample,  the  pump  beam  is  modulated  by  an  electro-optic 
modulator  (EOM)  at  a  frequency  from  a  few  kHz  to  a  few  MHz.  The  pump  beam  is  then 
absorbed  by  the  transducer  layer  (a  metal  thin  film,  usually  100  nm  Al)  deposited  on  the 
substrate,  creating  a  periodic  heating.  The  transducer  serves  as  the  sensor  at  the  same  time, 
because  the  reflectance  change  is  proportional  to  the  surface  temperature  rise.  To  monitor  the 
surface  temperature  change,  we  use  another  delayed  probe  beam.  The  probe  beam  passes  through 
a  mechanical  delay  stage  such  that  the  temperature  responses  are  detected  with  a  delay  time 
(<  8  ns  in  our  lab)  after  the  sample  is  heated  by  the  pump  pulse  using  a  photodetector.  The  signal 
from  thermoreflectance  change  is  then  extracted  by  the  lock-in  amplifier.  To  separate  the  probe 
beam  which  contains  the  reflectance  change  signal  from  the  pump  beam,  we  vertically  separate 
the  pump  and  the  probe  beam,  and  use  an  iris  to  block  the  reflected  pump  beam.  Figure  7  shows 
an  example  of  TDTR  signal  obtained  on  an  Al/GaN  sample.  The  thermal  transport  properties  can 
be  then  extracted  by  fitting  the  signal  with  the  heat  conduction  model. 
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Figure  6:  Schematic  of  TDTR  Measurement 


Figure  7:  Measured  TDTR  Temperature  Response  Ratio  (circles)  -Fin/Fout  of  bulk  GaN  at 
Pump  Modulation  Frequencies  of  2.03  MHz  and  5.03  MHz,  Coated  with  101-nm-thick  A1 

Transducer 

The  color  lines  are  predicted  by  a  theoretical  heat  transfer  model  with  thermal  conductivity 
kGaN  =  190W/mK,  CcaN  =  3.12  MJ/m3K,  and  interface  conductance  of  Gai-GuN  =  80MW/m2K. 

Figure  8  shows  our  calibration  measurements  of  the  TDTR  system  on  several  standard  samples 
including  diamond,  SiC,  GaN,  Si,  and  sapphire  with  thermal  conductivity  ranging  from  1  W/mK 
to  103  W/mK.  Our  measurement  agrees  well  with  the  reference  values. 
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fo  (MHz) 

Figure  8:  Calibration  Measurement  of  our  TDTR  System 

The  dash  lines  are  the  standard  reference.32  34 


4.2  Frequency  Filtering  Effect  of  the  Metal  Thin  Film  in  TDTR 

In  this  part,  we  discuss  the  physical  picture  of  the  frequency  filtering  effect  in  TDTR.  There  are 
two  major  factors  resulting  in  our  prediction:  i)  multi-frequency  heating  in  TDTR  and  ii)  finite 
response  time  of  the  metal  transducer  layer. 


Multi-frequency  Heating  in  TDTR.  In  TDTR,  the  femtosecond  pump  laser  pulses  are  modulated 
with  a  sinusoidal  envelope  (see  Figure  9a).  By  applying  Fourier  transform,  there  exist  a  series  of 
frequencies  coo  ±  ncos,  where  too  is  the  modulation  frequency  (the  sinusoidal  envelop  of  the  pump 
beam  in  Figure  9b).  When  the  aluminum  transducer  is  subject  to  the  multi-frequency  heating,  the 
transducer  would  respond  differently  to  the  heating  components  with  different  frequency.  This  is 
due  to  the  fact  that  there  exist  a  characteristic  response  time  in  the  Al/Substrate  system. 


(b) 


#  *  In-phiEi* 
O  D  O.rl-rrf-ph; 


F  1 


*  6  *  *  4  /(MHz) 


Figure  9:  (a)  Modulated  Pump  Pulses  in  Time  Domain  and  (b)  Pump  Pulses  in  Frequency 

Domain 
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Response  Time  of  the  Transducer/Substrate  System.  Due  to  the  existence  of  the  aluminum 
transducer,  it  takes  a  certain  amount  of  time  Tc  until  the  substrate  can  follow  the  temperature 
oscillation  at  a  certain  frequency.  When  the  periodic  laser  heating  incident  on  the  sample  surface, 
the  characteristic  time  Ti  it  takes  until  thermal  equilibrium  is  reached  in  the  aluminum  transducer 
can  be  approximated  by  equating  the  penetration  length  i  —  to  the  aluminum  thickness 

dAl,  P  1 


CAl  A2 
=i~dAi 

Km 


(16) 


where  Rai  and  Cm  is  the  thermal  conductivity  and  volumetric  heat  capacity  of  aluminum 
transducer.  After  the  transducer  film  reaches  local  equilibrium,  the  heat  conduction  is  dominated 
by  the  interfacial  heat  conduction,  and  the  temperature  oscillation  of  the  aluminum  film  can  be 
described  by  a  lumped-capacity  model: 


K\ldAl 


dZ 


At 


dt 


=  =  0)  -  r*) 


(17) 


where  G  is  the  interface  conductance  between  the  Tsubiz  =  0)  is  the  temperature  at  the  interface 
on  the  substrate  side.  The  characteristic  time  Ti  can  be  obtained  by  non-dimensionalize  equation 
(2): 


T2  —  C AldAl /G 

(18) 


The  characteristic  time  T2  describes  the  timescale  it  takes  until  the  substrate  can  follow  the 
temperature  oscillation  of  the  aluminum  transducer.  Therefore,  Ti  and  T2  combined  is  the 
characteristic  time  it  takes  until  the  temperature  of  the  substrate  is  following  the  periodic  laser 
heating. 


Tc  =  Tl  +  T2  = 


Mi 


(19) 


The  characteristic  time  rf  describes  the  timescale  it  takes  until  the  substrate  can  follow  the 
temperature  oscillation  of  the  aluminum  transducer.  Therefore,  Ti  and  T2  combined  is  the 
characteristic  time  it  takes  until  the  temperature  of  the  substrate  is  following  the  periodic  laser 
heating. 

Due  to  the  characteristic  time  introduced  by  the  transducer  layer,  the  substrate  can  only  follow 
the  temperature  oscillation  with  frequency  (x>h«  2u/tc,  and  the  temperature  oscillation  with 
frequency  higher  than  27T/Tc  would  be  filtered  by  the  transducer  layer  and  the  interface  such  that 
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the  substrate  would  not  respond.  The  physical  picture  is  shown  in  Figure  10.  Equation  (19) 
demonstrates  that  we  can  selectively  control  the  frequency  range  that  can  penetrate  into  the 
substrate  by  varying  the  thickness  of  the  transducer  layer.  The  thinner  the  transducer  layer  is 
deposited,  the  more  high-frequency  components  in  TDTR  experiment  would  be  allowed  to 
penetrate  into  the  substrate,  and  the  lower  the  thermal  conductivity  would  be  measured  due  to  the 
ballistic  effects  from  high  frequency  components. 

In  conclusion,  we  expect  that  the  measured  apparent  thermal  conductivity  would  be  a  function  of 
both  aluminum  thickness  and  the  modulation  frequency  /a(o>o,  cIai),  and  k.A  would  decrease  with 
increasing  modulation  frequency  coo  and  decreasing  cLai. 


Thick  Transducer,  high  Frequency 
responses  are  filtered  by  the  transducer. 


Thick  Transducer,  substrate  can  only  Thin  Transducer,  tow  frequency  response 
feel  low  frequency  heating 


Figure  10:  Physical  Picture  of  Frequency  Filtering  Effect  of  the  Transducer  Layer  at 

different  Heating  Frequencies 

4.3  Experimental  Observation  of  the  Frequency  Filtering  Effect 

Temperature-Dependent  Thermal  Conductivity  of  GaN.  Figure  1 1  shows  the  measured  thermal 
conductivity  of  GaN  from  80  K  to  300  K.  The  measured  thermal  conductivity  decrease  when  the 
modulation  frequency  increases.  This  might  be  attributed  to  the  quasi-ballistic  transport  when  the 
penetration  depth  dp  becomes  comparable  with  the  phonon  MFP  in  GaN.  We  treat  these  results 
with  caution,  because  using  an  erroneous  spot  radius  might  also  result  in  the  frequency 
dependent  thermal  conductivity.  To  rule  out  this  error,  we  performed  similar  thermal 
conductivity  measurement  on  sapphire,  which  is  more  sensitive  to  the  laser  spot  radius  in  TDTR 
experiment.  As  expected,  the  thermal  conductivity  of  sapphire  (35  W/mK)  is  frequency 
independent  at  room  temperature,  because  of  the  much  shorter  MFP  than  GaN.  However,  the 
frequency  dependence  starts  to  emerge  at  cryogenic  temperature.  For  example,  the  measured 
thermal  conductivity  of  sapphire  decreased  25%  at  100  K  when  the  modulation  frequency 
increased  from  1  MHz  to  10  MHz.  At  lower  temperature,  the  phonon-phonon  scattering  becomes 
weaker,  thus  the  MFP  becomes  longer  than  that  at  room  temperature.  Therefore,  the  phonon 
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ballistic  transport  and  frequency  dependent  thermal  conductivity  becomes  more  pronounced. 


(a) 


(b) 


Modulation  Frequency  (MHz) 


Figure  11:  (a)  Temperature  Dependent  Thermal  Conductivity  of  GaN  and  (b)  Frequency 
and  Temperature  Dependent  Thermal  Conductivity  of  Sapphire 

The  thermal  conductivity  is  observed  to  decrease  when  higher  modulation  frequency  is  used. 

Frequency-Al  Thickness  Dependent  Thermal  Conductivity  of  GaN.  Figure  12a  shows  the 
measured  thermal  conductivity  of  GaN  as  a  function  of  both  modulation  frequency  coo  and  the 
transducer  thickness  dAi.  Even  at  room  temperature,  the  ballistic  effect  is  pronounced.  For  each 
dAi,  the  thermal  conductivity  of  GaN  decreased  when  the  modulation  frequency  increases.  At  the 
same  time,  the  thermal  conductivity  also  decreases  if  thinner  aluminum  is  deposited  on  top  of  the 
GaN  substrate.  As  discussed,  the  thinner  aluminum  transducer  would  have  a  shortened  response 
time,  thereby  allowing  more  high  frequency  heating  components  penetrating  into  the  GaN 
substrate.  As  a  result,  the  high  frequency  heating  components  would  suppress  the  phonon  MFPs 
and  cause  a  lower  measured  thermal  conductivity.  We  have  also  designed  a  control  group.  We 
measured  the  thermal  conductivity  of  a  series  of  silica  samples  with  different  A1  thickness.  Since 
silica  has  very  low  thermal  conductivity  (~  1.4  W/mK)  and  short  MFPs  (tens  of  nanometers),  we 
should  not  observe  any  ballistic  effect  at  room  temperature,  and  the  thermal  conductivity  should 
be  independent  of  &>o  and  dAi. 
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Figure  12:  (a)  Temperature  Dependent  Thermal  Conductivity  of  GaN  and  (b)  Frequency 
and  Temperature  Dependent  Thermal  Conductivity  of  Sapphire 

The  thermal  conductivity  is  observed  to  decrease  when  higher  modulation  frequency  is  used. 


We  have  also  measured  the  frequency  and  A1  thickness  dependent  thermal  conductivity  of  GaN 
at  cryogenic  temperature  100  K.  Interestingly,  the  aluminum  thickness  dependence  is  even 
weaker  than  room  temperature.  Thermal  conductivity  is  decreased  by  only  10  %  when  dAi  is 
reduced  from  100  nm  to  20  nm,  as  compared  to  ~  20%  thermal  conductivity  reduction  at  room 
temperature.  The  reason  lies  in  the  temperature  dependent  thermal  conductivity  kAi  and  G  of  the 
Al/GaN  interface.  At  cryogenic  temperature,  both  k.Ai  and  G  decreased,  resulting  in  a  slower 
response  time  Tc=  Ca^m  +  Cai^ai  As  a  result,  it  becomes  even  harder  for  the  high  frequency 
heating  components  to  penetrate  into  the  substrate  due  to  the  stronger  frequency  filtering  effect. 


Figure  13:  Aluminum  Thickness  Dependent  and  Frequency  Dependent  Thermal 

Conductivity  of  GaN  at  100  K 
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Figure  14:  Temperature  Dependent  (a)  Thermal  Conductivity  of  A1  and  (b)  Conductance 

of  the  Al/GaN  Interface 
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5.  SUMMARY 


We  have  performed  first-principles  calculations  to  predict  the  phonon  and  thermal  properties 
The  phonons  with  mean  free  paths  larger  than  1  pm  contribute  ~50%  thermal  conductivity  for 
diamond,  GaN  and  SiC.  We  have  studied  dependence  of  GaN  thermal  conductivity  the  heater 
thickness  and  heating  frequency.  We  found  that  thick  heater  would  prevent  high  frequency  heating 
from  penetrating  into  the  substrate.  This  research  suggests  that  the  size  effect  of  thermal 
conductivity  in  GaN  devices  is  collectively  determined  by  several  characteristic  length  scales,  and 
the  previous  knowledge  of  size  effect  is  too  simple  to  understand  heat  dissipation  in  practical 
HEMT  devices. 
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LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


ACRONYM 

DESCRIPTION 

BTE 

Boltzman  transport  equation 

CVD 

chemical-vapor-deposited 

DARPA 

Defense  Advanced  Research  Projects  Agency 

DFPT 

density-functional  perturbation  theory 

DFT 

density  functional  theory 

EOM 

electro-optic  modulator 

GaN 

gallium  nitride 

HEMT 

high  electron  mobility  transistor 

MFP 

mean  free  path 

MTO 

Microsystems  Technology  Office 

PBTE 

Peierls-Boltzmann  Transport  Equation 

SiC 

silicon  carbide 

TDTR 

time-domain  thermoreflectance 
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